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In this note, we extend our previous work [Phys. Lett. B 644, 7 (2007), astro-ph/0609597 , and 
compare eleven interacting dark energy models with different couplings to the observational H(z) 
i 1 data. However, none of these models is better than the simplest ACDM model. This implies that 

^ | either more exotic couplings are needed in the cosmological models with interaction between dark en- 

ergy and dust matter, or there is no interaction at all. We consider that this result is disadvantageous 
C/2 to the interacting dark energy models studied extensively in the literature. 
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1 Introduction 

Nowadays, dark energy study has been one of the most active fields in modern cosmology CQ, since the 
discovery of the present accelerated expansion of our universe [21 02 HI [6l [7] . In the past years, many 
cosmological models are proposed to interpret this phenomenon. One of the important tasks is to confront 
them with observational data. The most frequent method to constrain the model parameters is fitting 
them to the luminosity distance 

d L (z)=(i+z) r 4k, 



/o H(z) 

which is an integral of Hubble parameter H = a/ a, where a = (1 + z) -1 is the scale factor (z is the 
rcdshift); a dot denotes the derivative with respect to cosmic time t. However, the integral cannot take 
the fine structure of H (z) into consideration and then lose some important information compiled in it (this 
point is also noticed in e.g. [S]). Therefore, it is more rewarding to investigate the observational H{z) 
data directly. 

The observational H(z) data we used here are based on differential ages of the oldest galaxies [9]. 
In [10j . Jimenez et al. obtained an independent estimate for the Hubble constant by the method devel- 
oped in [jj] , and used it to constrain the equation-of-state parameter (EoS) of dark energy. The Hubble 
parameter depends on the differential age as a function of redshift z in the form 

Therefore, a determination of dz/dt directly measures H (z) [11] . By using the differential ages of passively 
evolving galaxies determined from the Gemini Deep Deep Survey (GDDS) [H] and archival data [15] , 
Simon et al. determined H(z) in the range 0<z<1.8 CCD.- The observational H(z) data from [TT] are 
given in Table CO and shown in Figs.GHU 
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Table 1: The observational H(z) data CG31Q3] (see CLUE] also). 

These observational H(z) data have been used to constrain the dark energy potential and its redshift 
dependence by Simon et al. in [11] . Yi and Zhang used them to constrain the parameters of holographic 
dark energy model in CO)]. In [TI]> Samushia and Ratra have used these observational H(z) data to 
constrain the parameters of ACDM, XCDM and </>CDM models. Some relevant works also include C03 
OH QIl [HI [HI [20] for examples. 

By looking carefully on the observational H (z) data given in Table CO and shown in Figs. C3 C3 we 
notice that two data points near z ~ 1.5 and 0.3 are very special. They deviate from the main trend 
and dip sharply, especially the one near z <~ 1.5; the H(z) decreases and then increases around them. 
This hints that the effective EoS crossed —1 there. In our previous work [T5], we have confronted ten 
cosmological models with observational H (z) data, and found that the best models have an oscillating 
feature for both H{z) and effective EoS, with the effective EoS crossing —1 around redshift z ~ 1.5, while 
other non-oscillating dark energy models (e.g. ACDM, XCDM, vector- like dark energy etc.) cannot catch 
the main feature of the observational H (z) data. 

In Fig.CQ we show the quantity L(z) = H 2 (z)/Hq — il m o(l + z) 3 versus redshift z, which is associated 
with the fractional energy density of dark energy, for the fiducial parameters Hq = 72 km s _1 Mpc -1 and 
Qmo = 0.3 or 0.28, where the subscript "0" indicates the present value of the corresponding quantity. It is 
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easy to see that the fractional energy density of dark energy of the point near z ~ 1.5 is negative (beyond 
la significance). To avoid this, one can decrease the corresponding f2 m o or make the matter decrease with 
the expansion of our universe slower than a~ 3 . Inspired by this, it is natural to consider the possibility 
of exchanging energy between dark energy and dust matter through interaction. In fact, we considered 
the cases with constant coupling coefficient in |15j . However, we found that it is not preferred by the 
observational H(z) data. In the present work, we will explore more forms of couplings between dark 
energy and dust matter, in an attempt to find the couplings which can best describe the observational 
H{z) data. 



H =72 km/s/Mpc, n ra0 =0.30 



H =72 km/s/Mpc, n m0 =0.28 
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Figure 1: The quantity L(z) = H 2 (z)/H 2 - ft m0 (l + zf versus redshift z, for the fiducial parameters 
Hq = 72 km s _1 Mpc -1 and fl m o = 0.3 (left panel) or 0.28 (right panel). 



As extensively considered in the literature (see e.g. [2Il|22l2iaiH|2I3l2IIll22ll3l2SHffl 
[36l[39lllQlllIl|42l[43]), we assume that dark energy and dust matter exchange energy through interaction 
according to 

px + 3H {px + Px) = -SQHpm, (3) 
p ro + 3Hp m = 3QHp m , (4) 

which preserves the total energy conservation equation p to t + 3 if (ptot + Ptot) = 0. We assume that 
the EoS of dark energy wx = Px/px is constant, and consider a spatially flat Friedmann-Robertson- 
Walker (FRW) universe throughout this work. Notice that the coupling coefficient Q = Q(z) can be any 
function of redshift z. So, the interaction term 3QHp rn is a general form, in contrast to the first glance. 
Integrating Eq. ((4|), it is easy to get 

(5) 

where N = In a = — ln(l + z) is the so-called e-folding time; the constant proportional coefficient can 
be determined by requiring p m (N = 0) = p m o- Then, px can be also obtained by substituting p m into 
Eq. (|3|). From the Friedmann equation 3H 2 = &TrG(p m + px), the Hubble parameter is in hand. 

In the following sections, we will compare the observational H(z) data with some cosmological models 
with different couplings. We adopt the prior Hq — 72 km s _1 Mpc -1 , which is exactly the median value 
of the result from the Hubble Space Telescope (HST) key project [SI], and is also well consistent with the 



p m oc exp 



3(Q - l)dN 
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one from WMAP 3-year result [I]. Since there are only 9 observational H(z) data points and their errors 
are fairly large, they cannot severely constrain model parameters alone. We perform a x 2 analysis and 
compare the cosmological models to find out the one which catches the main features of the observational 
H(z) data. We determine the best-fit values for the model parameters by minimizing 

9 2 

X 2 {parameters) - ^ [Hmodjparameters; Zj) - H obs (zj)\ ^ 

where H mo d is the predicted value for the Hubble parameter in the assumed model, H a b s is the observed 
value, a is the corresponding la uncertainty, and the summation is over the 9 observational H{z) data 
points at redshift Zi. 



2 Simple couplings 

We firstly consider the simplest case, Q = 0, namely there is no interaction between dark energy and dust 
matter. It is easy to find that the Hubble parameter is given by [15] 



H{z) = H ^Jn m0 (l + zf + (1 - fi m0 )(l + z)3d+-x) > ( 7 ) 

where f2 m o = 8irG p m o / (3Hq) is the present fractional energy density of dust matter. By minimizing the 
corresponding x 2 , wc find that the best-fit values for model parameters are £l m o — 0.28 and wx = —0.90, 
while Xmin = 9-02 for 7 degrees of freedom and P (x 2 > Xmin) = 0-25. 

Next, we consider the case of Q — const.. By using Eqs. ([5|) and we can obtain the Hubble 
parameter as [15] 



By minimizing the corresponding x 2 , we find that the best-fit values for model parameters are fl m o = 0.72, 
w x = -3.70 and Q = 0.30, while Xmin = 8 - 48 for 6 degrees of freedom and P (x 2 > Xmin) = °- 21 - 
Although the Xmin 1S reduced slightly compare to the case of Q = 0, statistically this fitting has no 
significant improvement over the case of Q = 0, since the number of parameters is increased from 2 to 3. 
Another Q with only one parameter comes from the assumption 33, 35] I36j 



^ = ^V, (9) 

Pm PmO 

where £ is a constant parameter, which quantifies the severity of the coincidence problem. Following |33[ 
[35] [36], it is easy to find that the corresponding Q is given by (which is labeled as SCL) 

GO*) = i o +o° in k > ( 10 ) 

1 - S2 m0 + S2mo(l + Z)i 

where Qq = — (1 — fi m o)(£ + 3wx)/3; and the Hubble parameter reads 

h(z) = h (i + zf' 2 [n m0 + (i - n rnQ )(i + z)- £ ] - 3wx/{2 ^ . (ii) 

By minimizing the corresponding x 2 , we find that the best-fit values for model parameters are fl m o = 0.92, 
w x = —8.72 and £ = 1.49, while Xmin = 8 - 88 f° r 6 degrees of freedom and P (x 2 > Xmin) = 0-18. It is 
worth noting that the best-fit value Q m o = 0.92 is inconsistent with the results from clusters of galaxies [38] 
and 3-year WMAP g] etc. 
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We present the observational H(z) data with error bars, and the theoretical lines for these simple 
models with the corresponding best-fit parameters in Fig. [2] In fact, the cases of Q = and SCL cannot 
be distinguished significantly. It can be seen clearly from Fig. [3 that none of them may reproduce the 
sharp dip around z ~ 1.5; the data point near z ~ 1.5 deviates from model fitting by about 2a. In fact, 
it can be seen from Tabled all of these models are statistically even worse than the ACDM model. 




Figure 2: The observational H(z) data with error bars, and the theoretical lines for three simple 
models with the corresponding best-fit parameters. Left panel: the cases of Q = (solid line) and 
Q = const, (dashed line); Right panel: the case of SCL. 



3 More complicated couplings 

In this section, we will consider some more complicated couplings. Unlike in Section for these compli- 
cated cases, it is difficult to obtain the analytical px{z) from Eq. {3| with Eq. and then the Hubble 
parameter H(z). Therefore, the numerical methods [37] are necessary. For the convenience of numerical 
computing, we introduce p m = 87rGp m /(3i?o) an d Px = 8ttGpx /(3H^ ). Eqs. © and ([5]) become 

"" X - -3p x (l + w x )-3Qp m , (12) 



dN 



p m oc exp 



3(Q-l)dN 



(13) 



where the constant proportional coefficient in Eq. [|13p is determined by requiring p m (N = 0) = fl m ol the 
initial condition for integrating Eq. (fT2|) is px(N = 0) = 1 — fi m o- The Hubble parameter is given by 
H = Ha{px +pm) 1 ' 2 - ln the following subsections, we will consider some cases with parameterized Q(N). 
We restrict ourselves to the cases with only two parameters, since the data points are so few. To find the 
minimal x 2 efficiently, we scan the parameters space with a relatively large grid size at first, and then 
narrow the parameters space and scan with a small grid size. In the procedure, we impose the condition 
Px > 0, while p m > has been guaranteed by Eq. (fT3|) for real Q(N). 
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3.1 Linear coupling 

Here, we consider a linear coupling (LIN) 

Q(N) = Q + QiN, 
where Qo and Qi are constants. The corresponding p m is given by 



Q m0 exp 



3(Q - 1)N + -Q^ 2 



We hnd that the best-fit parameters are fl m o = 0.001, Wx = —0.80, Qo 
Xmin = 7 - 08 for 5 degrees of freedom and P (x 2 > Xmin) = °- 21 - 



-7.78 and Q 1 



(14) 
(15) 

-14.32, while 




Figure 3: The observational H(z) data with error bars, and the theoretical lines for the cases of LIN, 
COS1, COS2 and COS3, with the corresponding best- fit parameters. 



3.2 Couplings of cosine type 

We consider a coupling of cosine type (COS1) as 

Q(N) = A lC os(A 2 N), 



(16) 
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where A\ and A 2 are constants. The corresponding p m is given by 



Pm = ^raO exp 



-3^ + 3^ sin (A 2 N) 
A 2 



(17) 



We hnd that the best-fit parameters are f2 m o = 0.001, wx = —0.76, A\ = —6.09 and A2 = 2.87, while 



X/; 



7.07 for 5 degrees of freedom and P (x 2 > X 



2 ) 

rn i n I 



0.22. 



The next one (COS2) is given by 



It is easy to find that 



Pr> 



Q(N) = A\ cos (N/A 2 + Am) 



n mQ exp {-3N + 3AfA 2 [sin (N/A 2 + A x tt) - sin (^itt)]} 



-3.83 and A 2 



We find that the best-fit parameters are £l m o = 0.07, wx — —1.33, A\ 
Xmin = 5-89 for 5 degrees of freedom and P (x 2 > X m in) = 0.32. 
We consider the third case (COS3) 

Q(N) = Al cos (N/A 2 + ir/Ai) . 

The corresponding p m reads 

pm = to exp {-37V + 3A 3 X A 2 [sin {N/A 2 + n/Ai) - sin (tt/Ai)]} 



(18) 
(19) 

-0.04, while 

(20) 
(21) 



The best-fit parameters are found to be fl m o = 0.07, wx = —1.35, A\ = —4.41 and A 2 — 0.04, while 



X/, 



5.88 for 5 degrees of freedom and P (x 2 > X 



2 ) 

in in I 



0.32. 



In Fig. [21 we present the observational H(z) data with error bars, and the theoretical lines for the 
cases of LIN, COS1, COS2 and COS3, with the corresponding best-fit parameters. It can be seen clearly 
from Fig. [3] that the cases of LIN and COS1 cannot reproduce the sharp dip around z ~ 1.5. Although 
the cases of COS3 and COS4 have the oscillating feature, the data points near z ~ 1.5 and 1.75 deviate 
from model fitting beyond la. We notice that the best-fit parameter f2 m o for the cases of LIN, COS1, 
COS2 and COS3 are too small to be consistent with the results from clusters of galaxies [3H] and 3-year 
WMAP [4] etc. After all, even purely from statistical point of view, all these models are not preferred 
over the ACDM model. 



3.3 Couplings of Gaussian distribution type 

In this subsection, we consider some couplings of Gaussian distribution type. In fact, they can efficiently 
mimic the S function, namely, these Q can be rather large in a very narrow range of N and are approxi- 
mately zero for other N. 

We consider the case of original Gaussian distribution (GDI) at first, 



Q(N) = 



1 



V5FQ. 



■ exp 



(n - N s y< 



(22) 



(23) 



where Q s and N s are constants. The corresponding p m is given by 

[ 3 {N-N s 
p m = S2m0 exp -3N + - Erf I — — 

2 f x _ 2 

where Erf (x) = —= / e * dt is the well-known error function [37] . We find that the best-fit parameters 
Jo 



-Erf 
2 



are fL 



0.43, w x = -1.64, Q s = 0.007 and N s 



-0.85, while x 



2 

min 



5.89 for 5 degrees of freedom 



and P x 2 > XLJ = 0-32. 
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Then, wc consider a variant (GD2) of the original Gaussian distribution as 



(n - Ay 2 



and find that 



Erf 



<2s 

N — N s 



We obtain the best-fit parameters as Q m0 = 0.43, = —1.61, Q s — 0.014 and N s 
Xmin = 5-93 for 5 degrees of freedom and P (x 2 > Xmin) = 0.31. 



(24) 
(25) 

-0.89, while 
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Figure 4: The observational H(z) data with error bars, and the theoretical lines for the cases of GDI, 
GD2, GD3 and GD4, with the corresponding best-fit parameters. 



Also, other variant (GD3) of the original Gaussian distribution can be the form of 



Q(N) = - T = 7 - cxp 



(N-N a y 



Qt 



Its corresponding p m is given by 

p m = O to0 exp j-3Af + 3Q S 



)]}■ 



(26) 



(27) 
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We find that the best-fit parameters are £l m o = 0.47, wx — —2.10, Q s = 0.21 and N s = —0.83, while 
6.07 for 5 degrees of freedom and P (x 2 > Xmin) = 0.30. 
We consider the last one (GD4) as 



Xmin 



Q(N) 



(N-N s f 



It is easy to find that 



n m0 exp{-3N + 3Q 2 S 



■ exp 



We obtain the best-fit parameters as fi m o 

,i 



0.48 



0.42 and N s 



(28) 



(29) 



-0.82, while 



w x = -2.20, Q s 
Xmin — 6-43 for 5 degrees of freedom and P (x 2 > Xmin) = 0-27. 

In Fig. 21 we present the observational H(z) data with error bars, and the theoretical lines for the 
cases of GDI, GD2, GD3 and GD4, with the corresponding best-fit parameters. It can be seen clearly 
from Fig.UJthat the data points near z ~ 1.5 and 1.75 deviate from model fitting about la. Again, none 
of these models is better than the ACDM model. 



4 Conclusion and discussions 

In this note, we extend our previous work [15] . and compare eleven interacting dark energy models with 
different couplings to the observational H(z) data. In Tabled] we summarize all eleven models considered 
in this work. In addition, we show the results for the simplest ACDM model from [T5] together. Obviously, 
although the x m in °f au interacting dark energy models are lower than the simplest ACDM model, their 
Xmin/ do f are higher, and their P (x 2 > Xmin) are l° wer correspondingly. This implies that either more 
exotic couplings are needed in the cosmological models with interaction between dark energy and dust 
matter, or there is no interaction at all. We consider that this result is disadvantageous to the interacting 
dark energy models studied extensively in the literature (see [?T1 f^ ]. f^H I. f^H l f^T j f^H| . [^0] . f^T l [H^ l 

ESI ESI ESI HOI SH S3 H3] for examples) . 

Although the simplest ACDM model looks better, however, it is not the best model which is preferred 
by the observational H(z) data. In fact, as shown in our previous work [15j . the observational H(z) data 
favors the models which have an oscillating feature for both H(z) and effective EoS, with the effective 
EoS crossing —1 around redshift z ~ 1.5. Since as shown in the present work the interacting dark energy 
models fail, other physical mechanisms are needed to produce an oscillating feature for H(z) which can 
also fit the observational H(z) data fairly well. We leave this to our future works. 

We stress that other current data, such as SNe la [H El [7] , CMB [3] and so on, are not inconsistent 
with the ACDM model, as shown in e.g. [33]- Even for the observational H[z) data, we notice that from 
Table [H the ACDM model has Xmin = 9-04 for eight degrees of freedom with 34% probability, which is 
not unacceptably low. Before the new and improved H(z) data are available, it is too early to say that 
the ACDM model can be ruled out. 

In fact, we can also consider the couplings with three parameters, such as Q(N) = A\ 003(^42^ + ^4371-), 
Q(N) = A s exp [-(N - A s ) 2 /Q 2 ] , Q(N) = A s exp [-(A - N s f/Qf\ and so on. However, they are less 
attractive, since the data points are so few. Also, in addition to the couplings of Gaussian distribution type, 
Q(N) can mimic S function through e.g. Q(N) = Q h scch[2(N - N h )], Q(N) = Q h sech[2Q h {N - N h )], 
Q(N) = Qh sech[2Q^(N — Nh)] and so on, where sech(x) = 1/ cosh(a;). But we do not consider these 
cases in this work any more, since it is expected that they are physically similar to the cases of Gaussian 
distribution type. 

If we discard the condition px > 0, the fit can be improved, such as Xmin = 4.63, Xmin/dof = 0.93 
and P (x 2 > Xmin) = °-46 for parameters O m0 = 0.02, w x = -1.73, A\ = 3.16 and A 2 = 0.03 for the case 
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1.21 
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Table 2: Summarizing all eleven models considered in this work. In addition, we show the results for 
the simplest ACDM model from [TH] together. We have adopted Hq = 72 km s^ 1 Mpc" 1 exactly in the 
computations. 

of COS2, while Xmin = ^.77, X m i n /dof = 1-15 and P (x 2 > Xmin) = u -33 for parameters Q m0 = 0.57, 
wx = —2.41, Ai — 1.35 and A^ = 0.09 for the case of COS3. However, we consider that they are not 
physically acceptable, since their px become negative for some intervals of redshift z. 

It is worth noting that for the cases with couplings of Gaussian distribution type, careful treatment 
in the numerical computing is necessary. If one naively uses the Bulirsch-Stoer method or the adaptive 
Runge-Kutta method [37] to integrate Eq. (fl"2"|) directly, misleading results can be obtained, mainly due 
to the (5 function feature of the Gaussian distribution. For instance, one may find that Xmin = 1-64, 
Xmm/dof = 0.33 and P (x 2 > Xmin) = O- 90 for parameters fl m0 = 0.41, w x = -1.42, Q s = 0.10 and 
N s = —0.89 for the case of GD3. However, this is misleading. We present the corresponding results of 
H(z) and px, p m in Fig. [5j Obviously, the unaccounted fiber- like dips are not physical results. By using 
a more delicate treatment in the numerical computing, the right results can be obtained, as presented in 
Section |U and Fig.H 
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x 2 
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^ (X ^* Xmin) 
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1.29 


0.25 


Q = 


8.88 


1.48 
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Q = const. 


8.48 


1.70 


0.13 


SCL 


8.81 


1.76 


0.12 



Table 3: Summarizing the four models in which Hq is considered as a free parameter with the prior 
H = 72 ± 8 km s" 1 Mpc" 1 at 1 a uncertainty [33]. 

After all, it is worth noting that we adopt the exact Hq = 72 km s _1 Mpc _1 throughout this work. 
However, we admit that the uncertainty in Ho should be taken into account. For instance, it is suggested 
that in [33] H = 72 ± 8 km s^Mpc" 1 at la uncertainty, in [35] H = 68 ± 7 km s^Mpc -1 at 2a 
uncertainty, and in 46J Ho = 62.3 ± 6.3 km s" 1 Mpc" 1 at 1 a uncertainty. So, we should examine the 
effects of the uncertainty in Ho on our conclusions. Here we consider the four simplest models studied in 
this work for examples. They are the models ACDM, Q — 0, Q — const, and SCL. The Ho is considered 
as a free parameter in these models, with the prior Hq = 72 ± 8 km s _1 Mpc -1 at 1 a uncertainty [33]. In 
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this case, the total x 2 should be the summation of the one in Eq. © and Xho = (-Ho — 72) 2 /8 2 . We find 
that the best-fit values for model parameters of ACDM are fl m o — 0.31 and Hq = 71.44 km s" 1 Mpc -1 , 
while Xmin = 9-03 for 7 degrees of freedom and P (x 2 > Xmin) = 0-25. For model Q = 0, the best-fit 
parameters are Vl m o = 0.18, wx — —0.45 and H = 68.53 km s _1 Mpc -1 , while Xmin — 8-88 f° r ® degrees 
of freedom and P (x 2 > Xmin) — 0.18. For model Q = const., the best-fit parameters are Q m o = 0.72, 
w x = -3.68, Q = 0.30 and H = 71.88 km s -1 Mpc -1 , while x mi „ = 8.48 for 5 degrees of freedom and 
P (x 2 > Xmin) = 0-13. For model SCL, the best-fit parameters are Q m0 = 0.94, w x = -8.45, £ = 0.88 
and H = 69.79 km s -1 Mpc -1 , while x mm — 8-81 for 5 degrees of freedom and P (x 2 > Xmin) = 0.12. 
We summarize these results in Table [3] Comparing with the corresponding results of these four models 
in Table [2] where Hq is taken as 72 km s _1 Mpc _1 exactly, we find that the uncertainty in Ho cannot 
significantly improve the fits. Thus, even the uncertainty in Hq is taken into consideration, our conclusions 
remain unchanged. 




Figure 5: The numerical results for the case of GD3 with the corresponding "best-fit" parameters, which 
come from a careless numerical computing. See text for details. 
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